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_ O ' In 1960 Edward Lorenz (1917-2008) published a pioneering work on the 'maximum simpli- 

. fication' of the barotropic vorticity equation. He derived a coupled three-mode system and 

interpreted it as the minimum core of large-scale fluid mechanics on a 'finite but unbounded' 
domain. The model was obtained in a heuristic way, without giving a rigorous justification 
for the chosen selection of modes. In this paper, it is shown that one can legitimate Lorenz' 
choice by using symmetry transformations of the spectral form of the vorticity equation. 
The Lorenz three-mode model arises as the final step in a hierarchy of models constructed 
'■^ . via the component reduction by means of symmetries. In this sense, the Lorenz model is 

^ ' indeed the 'maximum simplification' of the vorticity equation. 

q' 

O 

'c/3 ■ 1 Introduction 

Q_il Symmetry is one of the most important concepts in numerous branches of modern natural 

science. The exploitation of symmetries of dynamical systems may lead to a more efficient 
ff^ . treatment of the differential equations describing these systems via a reduction of the information 

I that is necessary in order to account for model dynamics. 

\^ I There are many ways to utilize symmetries of differential equations, including the systematic 

■ construction of exact solutions of PDEs, determination of conservation laws and construction 

, of mappings that relate or linearize differential equations (see e.g. [Il[8l[9]). Also, there is some 

I work on symmetries in the study of dynamical systems and bifurcation theory [3], giving rise 

' to the study of equivariant dynamical systems. In the present work, we will take a related but 

. somewhat different direction to investigate how the Lorenz-1960 model [7] can be derived in a 

rigorous way. 

In this classical work, Lorenz considered the spectral expansion of the barotropic vorticity 
equation on a torus. In what follows, he sought for the minimum system of coupled ordinary 
first order differential equations for the Fourier coefficients that is necessary to still account 
for the nonlinear interaction of modes. In doing so, he first restricted the range of indices in 
the infinite Fourier series by the values {—1,0,1}. His crucial step to achieve the maximum 
simplification was the observation that all but three of these remaining coefficients retain their 
particular values once they are taken. We aim to give a justification of this observation by 
interpreting it as a condition of symmetry. To be more precise, we show that this simplification 
is possible due to the corresponding spectral counterparts of point symmetry transformations 
of the vorticity equation in physical space. These symmetries are preserved under truncations 
of the infinite Fourier series and are inherited by the spectral set of equations for the Fourier 
coefficients. 

A justification of the selection of modes in finite-mode models by considering inherited sym- 
metries may be potentially useful also in more general situations. It could provide an additional 
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criterion addressing the important question which modes in the reduced model should be re- 
tained and which may be neglected. In this sense, using induced symmetries may supply special 
kinds of truncations that are designed to preserve e.g. the invariants of the parental model also 
in the truncated dynamics (e.g. [lO]). This may lead to more concise and consistent finite- mode 
representations. As stated above, this approach is to be distinguished from the field of equiv- 
ariant dynamics. The main difference is that in the latter usually no exhaustive and rigorous 
calculations of symmetry groups are given. 

The organization of this paper is the following: Section[2]is devoted to discrete and continuous 
symmetries of the barotropic vorticity equation in a non-rotating reference frame. We expand 
the vorticity in a double Fourier series and discuss how the symmetries of the equation in physical 
space are induced to symmetry transformations in terms of the Fourier coefficients. In section [3] 
the initial model of eight ODEs will be presented amongst a discussion of the induction of the 
symmetries to this truncated system. Finally, in section H] subgroups of the whole symmetry 
group of the truncated spectral vorticity equations will be used to derive hierarchies of reduced 
models describing the evolution of the relevant Fourier coefficients. The main result is that 
the Lorenz-1960 model is indeed the maximum simplification of the dynamic equations. It is 
the minimal system that can be obtained by using induced discrete symmetry subgroups of the 
spectral vorticity equation for component reduction. 



2 Symmetries of the barotropic vorticity equation 

The inviscid barotropic vorticity equation in an inertial system in stream function form reads 

|v^^ + 2*«vV-|*Av=* = o. (1) 

ot OX ay oy ox 

where ip is the stream function generating two-dimensional nondivergent flow in the (x, y)-plane. 
It states the individual or Lagrangian conservation of the vorticity 

C = vV- 

Eqn. ([!]) possesses the eight-element group of discrete symmetries, generated by the elements 

ei : {x, y, t, ip) (x, -y, t, -ip) 
62: {x,y,t,'ip) {-x,y,t, -ip) 
63 : {x, y, t, ip) {x, y, -t, -ip). 

Note that these transformations are involutive (i.e. = 1, i = 1,2,3) and, moreover, they 
commute (i.e. eiCj = ejCi, i,j = 1,2,3, i ^ j). Each 6j generates a copy of Z2, the cyclic group 
of order 2. The group of all discrete symmetry transformations of the vorticity equation may 
then be written as Z2 ® Z2. 

For sake of completeness, we also present the basis operators of the maximal Lie invariance 
algebra of ([1]), which were computed using the program LIE by A. Head [6]. They read 

^ dt ^ ^ dy ^ dx 2 ^ ^ ^ dip 

^ ^ dx ^ dy dip ^ dt dip^ 
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where /, g and h are arbitrary smooth functions of t. Thus, in addition to the discrete symmetries 
ei, i = 1,2,3, the vorticity equation ([T|) possesses the time translations (generated by v^), the 
rotations with constant velocities and on constant angles (v„ and v^, respectively), the gauging 
of stream function with arbitrary summands depending in t {Z{h)), translatory motions with 
arbitrary (nonconstant) velocities (<^i(/) and ^^2(5)) and (separate) scaling of the space and 
time variables (Pi and 'D2)- In contrast to the other basis operators, the operator v„ has no 
counterpart in the three-dimensional case and leads to nonlocal transformations in terms of the 
fluid velocity and the pressure. This singularity, from the symmetry point of view, of the two- 
dimensional vorticity equations in terms of the stream function was first observed by Berker [1] 
and later re-opened (c.f. [3]). In this paper we simultaneously use discrete symmetries and some 
of transformations from the continuous symmetry group. 

We now expand the vorticity in a double Fourier series on the torus, 

C = Cm exp(im ■ x), x = xi-|-yj, m = mii + m2j, m = miki + m2lj, 

m 

where k and I are nonzero constants, i = (1,0,0)"'", j = (0, IjO)""", mi and m2 run through the 
set of integers and the coefficient cqo vanishes. Inserting this expansion in the vorticity equation 
([I]) gives its spectral form [7]: 

^ = - L rii^2 (k- [m X m]) , (2) 

where k = (0,0, 1)""". The transformations ej induce transformations of the Fourier coefficients. 
Thus, in spectral terms the action ({x, y, t) = —C{x, —y, t) of ei on ( has the form 

Cm exp(im • x) = - ^ Cm exp(i(miA;j; - m2ly)) = - ^Cmi,-m2 exp(«(mi/cx -|- m2ly)), 

mm m 

upon changing the summation over m2. Hence as a consequence of 

the x-reflection ei. Finally, these and similar computations give the transformations 



ei : 
62: 



-"17117712 '~-77l\,—77l2 



-mim2 



C3 ■ Cmim2 ^ Cjmr?i2; ^ * ^• 

It is obvious that the induced transformations for the coefficients Cm are symmetries of sys- 
tem 

Some of continuous symmetries of the vorticity equation also induce well-defined symmetries 
of system ([2]) which are applicable within the framework of our approach. Thus, the rotation 
on the angle vr coincides with 6162. (In fact, only one of the reflections ei and 62 is independent 
up to continuous symmetries.) More nontrivial examples are given by the space translations, 
which induce the transformations 

n ■ c J^ike 

ye • '-mim2 ^ '-'mim2 

n ■ r p«"i2fe„ 

ye- '-mir?i2 '~'77%\77i2- 

For the values e = 7r//c and e = vr/Z, respectively we have ^7^/^ • Cmim2 ~^ (— l)'"^Cmim2 and 
<1-kIi'- Cmim2 ~^ (" l)™^Cm^m2. 1^°^^ ^he Sake of brevity, we will set p = p^/;. and q = q^ji. It is 
easy to see that the transformations p and q are involutive and commutes with ei, 62 and 63. 
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3 Discrete symmetries of the truncated system 



The restriction of the range of indices in ([2]) by {—1,0,1} x {—1,0,1} leads to the following 
eight-mode model related to the vorticity equation: 
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klc-iQCQi = \y-^2r—p - -j^j [c-ioci -1 - cioc_i „i] . (3) 

System ([3]) consists of first-order ordinary differential equations. Hence the problem on descrip- 
tion of its point symmetries is even more difficult than its complete integration (c.f. [8]). At 
the same time, some symmetries of ^ are in fact known since they are induced by symmetries 
of the vorticity equation. In particular, the symmetric truncation guarantees that the resulting 
system ([3]) inherits the discrete symmetries of the initial system ([2]), induced by ei, 62, 63 and 
their compositions. (This is an argument justifying such kind of truncation.) Any truncation 
also preserves the symmetries Pe and q^, in particular, p and q. 

Induced transformations are of crucial importance for deriving the Lorenz system. 



4 Component reduction of the truncated system 

The transformations ei, 62, p and q and their compositions act only on dependent variables 
and, therefore, can be used for component reductions of system ([3]) within our approach^ The 
technique applied is similar to that developed for invariant solutions without transversality (c.f. 
[21 [9]). However, it is not apparent that the element 63 can be used for this purpose. We thus 
restrict ourself to the group G ~ Z2 © Z2 © Z2 ffi Z2 generated by ei, 62, p and q. We will apply 
the following subgroups of G: 

5i = {l,ei}, 51 = {1,62}, 52 = {l,p}, S'^ = {l,q}, 
'S'3 = {l,pe2}, S'3 = {l,gei}, 5'4 = {1,m}, 
5*5 = {l,]5ei}, 55 = {1,^62}, 56 = {1,6162}, 

5*7 = {I,j56i62}, 57 = {1,^6162}, 5*8 = {l,pg6i}, 5*8 = {l,pg62}, 
59 = {I,j5g6i62}, 5*10 = {l,PQei,pge2,6i62}, 

Sii = {I,pei,qeie2,pqe2}, S'n = {1,^62,^6162,^61}, 

5l2 = {I,p6l,^62,]5g6l62}. 

By 1 we denote the identical transformation. There are yet other subgroups of G but as we will 
see they do not lead to nontrivial reduced systems. 

Note that the subgroup S[ will result in the same reduced system as 5*1 (up to the re-notation 
{x,k) ^ (y, 0)- ^ similar remark holds also for all the subgroup marked by prime, so we only 
have to consider reductions with respect to the subgroups without prime. This should be done 
subsequently. 

^See the appendix for a depiction of these transformations. 
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4.1 Trivial reductions 



The tuple (cjj) is invariant with respect to the transformation ei if and only if the following 
identifications hold: 

C-1,-1 = — c„ii, Co^_l = — Cqi, Ci^_i = — Cii, Cio = — Cio, c_io = — c„io. 

The last two conditions require that cio = c_io = 0. Inserting these identifications in ([3]) leads 
to a trivial system only. Hence, the subgroup 5i cannot be used for a component reduction 
of ([3]). As a consequence, no subgroup of G that contains the transformation ei can be used for 
this purpose. 

The same statement is true for the transformations 62, j5, q, pe2, qei and pq. Finally, any 
subgroup of G containing one of the elements ei, 62, p, q, pe2, qei or pq gives a trivial reduction. 



4.2 Reductions in three components 

The tuple (cij) is invariant under the subgroup Ss generated by the transformation pqei if and 
only if the following equalities hold: 

Co -1 = coi, c_i _i = -c_ii, ci _i = -cii. 

This transformation allows us to reduce ([3]) to the nontrivial five-component system 

dcio „ f 1 1 \ , , 



dt \k'^ + P P 

dc_io . / 1 1 



dt \k^ + P P 

dcoi / I 1 



dcii 




1 


dt 


(} 


" ¥ 


dc_ii 






dt 




~P 



klc^iicoi —n~ = ( 12 ~ 72 ) klc^iocoi (4) 

fc/ [cioc_ii - CiiC_io] . 



dt V't^ + ^^ 

The subgroup leads to a similar reduction under equating 

co,-i = -coi, C_i_i = C_ii, Ci „i = Cii 



4.3 Reductions in four components 

Utilizing the subgroup Sq = {1,6162} by means of similar consideration as in the previous 
sections leads to the following nontrivial reduced system of ([3]) after equating c_ii = ci^_i, 
c-1,-1 = cii, 60,-1 = 601 and e_io = eio: 

^^ciocoi ^ = ( ^2^^2 - ^) ^^[^11 - ci,-i]coi 

^ ) A;/ [ci _i - en] cio- (5) 
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Similar four-component reductions are also given by the non-primed subgroups 5*7 (e_ii = 
-ci-i, e_i _i = -en, cq-i = eoi and e_io = -eio) and Sg (e_ii = ci_i, c_i_i = en, 
60 -1 = -601 and e_io = -eio). 
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4.4 Maximal reduction: subgroup S^o 

Finally, let us derive the reduced system associated with the subgroup Siq . This reduced system 
coincides with the Lorenz-1960 model. Since all transformations mutually commute, there are 
three equivalent ways for the reduction. The first way is to perform the reduction in a single 
step from system ([3]) using the subgroup 5io. Because the system ([5]) possesses the symmetry 
transformation induced by pqei, we can also start with ([5|) and compute the induced component 
reduction. Alternatively, it is possible start with (j4]) and apply the transformation induced by 
6162. This will give the same reduced system associated with the subgroup 5*10 • 

The subgroups Sn and 5i2 lead to similar reductions connected with the Lorenz' reduction via 
transformations generated by p^/2fc (77^/2/- Therefore, up to symmetries of ([3]) induced by the 
symmetries of the vorticity equation, there is a unique nontrivial reduction in five components 
using subgroups of G. 

However, let us note before that we have already employed the transformation 

{x,y,t,'il;) {-x,-y,t,Tp) 

which accounts for the identification c-m = Cm- Observe that the connection between the 
complex and real Fourier coefficients is given by 

Om — Cm -|- C_ni) — ^(Cm C_m)' 

Hence we can already set 

_ 1 

Cm — '^dm- 

This explains the first observation by Lorenz : If the imaginary parts of the Fourier coefficients 
vanish initially, they will vanish for all times. Mathematically, this is justified since we have 
used the symmetry transformation group 5*6. 

Now employing the transformation pqei at this stage of simplification, we have the identifi- 
cation 

ci _i = -cii <^ ai _i = -an. 

This is the second observation by Lorenz [7], namely that if ai^_i = —an initially, then the equal- 
ity will hold for all times. That is, Lorenz heuristically discovered the symmetry transformation 
pqei to obtain his famous finite-mode model. 

It is straightforward to see that the resulting system leads to the Lorenz system upon setting 
A = aoi,F = aio, G = ai _i: 

dt k^ + Pj 

= (/2 - ■^^) ^^^^ 

dt - 2[p k^)''^^^- 

The Lorenz-1960 system is thus the result of invariance of the system ([3]) under the subgroup 
SiQ. Hence it is truly the maximum simplification that can be obtained within our symmetry 
framework. We have exhausted all the symmetry groups consisting of chosen discrete symmetry 
transformations and leading to nontrivial reduced systems. 



6 



Acknowledgments 



The authors are grateful to the referees for their remarks. The authors thank Michael Hantel 
and Leopold Haimberger for helpful discussions. The research of ROP was supported by the 
Austrian Science Fund (FWF), START-project Y237 and stand-alone project P20632. AB is a 
recipient of a DOC-fellowship of the Austrian Academy of Science. 

Appendix: Depiction of symmetry transformations 

For the readers' convenience, we graphically illustrate some of the symmetry transformations 
used for the component reduction. An arrow between two components means the identification 
of the respective coefficients. If there is a minus sign next to an arrow, the corresponding 
identification is up to minus. For any arrow starting and ending in the same coefficient, the 
identification up to minus means that this coefficient vanishes. The coefficients remained in the 
finite-mode models after the identification are displayed in bold. 
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